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In the framework of a transient accretion disk at the core of a gamma-ray burst we compute 
possible periods of Lense-Thirring precession. Next, we evaluate the putative gravitational waves 
associated with such dynamical setup. Assuming a characteristic time-profile for the gamma-ray 
emission of a disk-jet system, we obtain light-curves presenting a time microstructure similar 
to that reported in some GRB events. After adjustment of the parameters out of two specific 
GRBs we evaluate the detectability of the gravitational waves produced by the precession of this 
accretion disk. As a conclusion, our analysis shows that some GRBs are likely to be probed with 
Advanced LIGO. 
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Gravitational waves from precessing engines in GRBs 



1. Introduction 

It is widely thought that gamma-ray bursts (GRBs) are stellar events produced out of a central 
object that consists of a massive black hole with a transient, hot and dense accretion disk. Such an 
engine might result from a failed supernova or coUapsar (e.g. [|T]]) or the merging of two compact 
objects (e.g. [Q]). Apparently, the accretion of matter onto the black hole is the energy source of 
the system. Part of this energy is ejected through an ephemeral relativistic bipolar jet, normal to the 
midplane of the disk, radiating an extremely intense burst of gamma-rays. The duration of these 
gamma-ray pulses ranges from milliseconds up to minutes. This picture is expected to underly 
most of the observed GRB light-curves. 

The usual explanation for the temporal structure of such light curves relays on the formation 
of shock waves that convert bulk kinetic energy into relativistic energy of the particles in the jets. 
Charged particles then cool by synchrotron and inverse Compton emission. The shocks can be 
either internal to the jet and produced by colliding shells with different Lorentz factors (e.g. |[ 
^) or the result of interactions with the external medium (e.g. [^]). 

Nevertheless, some peculiar GRB light curves are difficult to explain by means of shocks only, 
particularly those having slow rises and fast decays (e.g. [^). In this respect, it has been suggested 
that the precession of the jet might play an important role in the formation of its time microstructure 
and this could be valid in both long and short gamma-ray bursts (e.g. 10, ll|]). In [[111] amodel 



for precessing jets has been based on spin-induced precession of a neutrino-cooled massive disk. 
We will use this model to obtain information about the inner system from the GRB light-curves. 

Our goal here is the study of another expected byproduct of disk precession, namely the pro- 
duction of gravitational waves. Strong gravitational radiation is likely in several situations: when 
the gravitational collapse originating the burst is non-spherical, in the presence of strong inhomo- 
geneities in the accretion disk, or, for short GRBs, as a result of the spiral merging of compact 



objects (e.g. [ |12| , |13||). 

In the next sections we will present a model for the production of a gamma-ray light curve by 
a central engine consisting of a massive black hole, a transient accretion disk and a jet. Considering 
Lense-Thirring precession we will then compute gravitational waves emitted by such a source and 
finally discuss their possible detection. 



2. Accretion disk and spin-induced precession 

In any model for a transient accretion disk formed after massive collapse or merging of com- 
pact objects, it is reasonable to consider it to be initially misaligned with the resulting black hole. 



As first noticed by Lense and Thirring in the general case [ |14| | , in the central engine of a GRB this 
misalignment will cause disk precession. This phenomenon is originated in the dragging suffered 
by the inertial frames near a rapidly spinning black hole. 

When the Mach number inside the disk is below 5, which seems to be the case in most GRBs, 
it is likely that the disk precesses approximately like a rigid body, i.e., it presents no warping [p3|]. 
Of course, the precession of the disk should lead to the precession of the jets, yielding a likely 



source of temporal micro- variability in the gamma-ray signal [11]. 
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The typical disk accretion rate, ranging from ~ 0.1 to lOMg, is expected to significantly 
vary in the outer part of the disk. In contrast, for the inner disk a constant accretion rate should 
be a valid approximation (e.g. [1^, 17]). The conservation of mass, falling with a velocity v,- ~ 



r^J GMh^r ^ at a distance r from the black hole axis, is given by 

M = -2nrv,.L{r) 



(2.1) 



where £(r) = 2p{r)H{r) is the surface density, H{r) is the disk half-thickness, and p(r) is the 
mass density of the disk. Now, conservation of angular momentum and energy can be used to 
numerically compute r(r) and H{r), considering that the heat generated by friction can be balanced 
by advection and neutrino emission [[III]. 

Without nutation, the precession period of the disk Tp can be related to its surface density by 



[18, 11] 



Zp = / — smt^a^ = 27rsmt^ — 
Jo 7d 7d 



where the values of the disk angular momentum and applied precessional torque are 

^.d = 271 / I(r)ak(r)r^ dr, 



Td = 47r^sin0 / r(r)ak(r)Vp,e(r)r' dr. 

•'Rms 



The relativistic Keplerian angular velocity reads 

.3 



c 



3/2 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



where a* is the spin parameter, Rg = GMbb/c^ the gravitational radius, and 



271 



1 



(2.6) 



is the nodal frequency resulting from the perturbation of a circular orbit in the Kerr metric [|^. 
Since the action of viscous torques leads to the alignment of the very inner part of the accretion 
disk with the black hole equator plane [^, the precessing part of the disk ends at an outer radius 
/?out> extending from an inner radius R^s = i^ms^g. where 



with 



and 



U = 3+A2T[(3-Al)(3+Al+2A2)]^/^ 
Ai = 1 + (1 -a2)l/3[(i +^^)l/3 + (1 _^^)i/3] 

A2 = {3al+Ajy^\ 



(2.7) 



(2.8) 



(2.9) 



The minus sign in i^ms corresponds to prograde motion (a* > 0), whereas the plus sign corresponds 
to retrograde motion (a* < 0). 
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3. Gravitational waves from Lense-Thirring precession 



The gamma-ray luminosity produced by a relativistic jet depends on the observation angle i//^ 
as given by [P] 



L{^if) 



27 



-0.6*(x) 



d'^ix) 
dx 



(3.1) 



Here, \\r is the angle between the jet and the observer, x = lOsin and ^{x) = ^ ln(l + Ax^). The 
intrinsic time dependence of the light signal, on the other hand, might be characterized by a FRED 
(Fast Rise and Exponential Decay) function. 



I{t) =Ni[l-e 



71 

tan 

2 



-1 



(3.2) 



where Nj is a normalization constant such that the maximum of the signal corresponds to unity. 
The initial rise, plateau, and decay timescales are Tnse, Tpiat, and Tjec respectively. 

In order to make our predictions we shall consider two specific GRBs, the short GRB 990720 
and the long GRB 990712. For these events, we use F{t) = I(t) x L{\\r{t)) to reproduce the ob- 
served light curves. Here, the angle between the jet and the observer is time-dependent because 
the azimuthal angle of the jet is (?) =2K{t/ Tprec ) as a result of the disk/jet precession (see Fig.p. 
Both angles are related by 



cos l/A = rjet • ^obs = cos sin a sin 0obs + cos a cos 0obs • 



(3.3) 



Taking this into account, we can obtain the relevant timescales and precession period by the proce- 
dure above described. In both cases, the observer is located at 0obs = 2° with respect to the z-axis 
perpendicular to the black hole equator, and we set 0obs at = (see Figjl]). 

3.1 Emission of gravitational waves 

Axis-symmetric bodies (i.e. with inertial moments I\ = h) in precession emit gravitational 



waves with an amplitude given by pl| , p^ ] 
where 



with 



hj^{t) = F+.i cosHf + F+^2Cos2nf, hy^{t) = Fx ^ sinHf + Fx,2sin2nf, 



F+i = /iq sin 2a sin I cos I, F+2 = 2/zosin^a(l +cos^i) 



(3.4) 



(3.5) 



x,l 



/iq sin 2a sin I, f x ,2 = 4/jo sin^ a cos i , 



and h'n 



{h —h)^/d. Here, a is the angle between the angular momentum of the disk and 



that of the black hole (see Fig.|I|), i is the angle between the z-axis of the detector and the signal 
direction of arrival (line of sight), and d is the distance to the radiating body. The principal moments 
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of inertia are 73 = (x^ +y^)p{r), I\ = /y(z^ +y^)p{r). The frequency of the gravitational waves 
(GW) are fi = Q./ {In) and /2 = 2Q./ {2k), which are related to the angular momentum of the body 
by n = L//i. 

Since we are dealing with bursting events, the GW signal is expected to be significant for a 
brief time Tpiat- We therefore modulate the signal given by Eq. ( |3.4[ ) with a gaussian 



h{t) = Vec(0 ^ 



' plat 



(3.6) 



as usually done to describe GW signals from bursting sources (e.g. [ ]22| , |23|]. The angular fre- 



quencies that contribute to the waveform of Eq. ( [3.6[ ) can be better obtained looking at the Fourier 
transform of Eq.(||^ h{(o) = h+{(o) + h^{(o), where 



h+((0) 



Vat 




and 



hxico) 



IT, 



plat 



x,l 



e 2 



+ e~ 



((B-£l) 



,2 e 



2 



(£0+2^)2 



-^(co-2a)2 



From these expressions it is apparent that the main contributions to the signals are close around Q. 
and 2Q., and the width of the interval is Tpj^^. Given the typical durations of GRBs, the frequency 
spread is narrow for all burst except for those with durations much shorter than 1 s. 
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Figure 2: (up) GRB 990720 light curve, experimental (left panel) and model obtained (right panel), (down) 
GRB 990712 light curve, experimental (left panel) and model obtained (right panel). 



3.2 Numerical predictions of Lense-Thirring gravitational waves 



In order to judge detectability, we compute the root-sum-square amplitude (e.g. [22, 23 1) using 



Eqs. (3.5) 



hUf) = \ dt{h\+hl)e 



plat 



(3.7) 



For illustration, we choose i = 45° and a = 10° , 20°. In Fig. |3| we plot the root-sum-square ampli- 
tude /irss, as a function of one of the resonant frequencies, f = Q./{2n), and also as a function of the 
outer radius of the inner precessing disk. In the figure we include the expected sensitivity for Ad- 
vanced LIGO [p4|]. The parameters used for finding r(r) and H{r) are Mbh = 3M0, a* = 0. 1, a vis- 
cosity parameter a = 0. 1 , and the different mass loss rates are M = {0. \Mq s~ \ IM0 s^ ^ , lOM© s^ ' }. 
The distance to the GRB is taken as = 100 Mpc. 

Fig. § makes clear that there are better chances of detection for accretion rates higher than 
IMq s^' and outer radii between 10^ and 10^ cm. Note that when the accretion rate is very high 



the disk may become advection dominated rather than cooled by neutrino emission [ ]25| ] but this 
should not affect the dynamics in the gravitational field. On the other hand, large accretion rates 
can be only supported in long GRBs, so we conclude that there is only a good prospect for detection 
of gravitational waves from precessing disks of nearby {d < 100 Mpc) long events. Such events 
are likely related to the death of very massive stars, so the host galaxies should have active star 
forming regions. 

Two low-luminosity long GRBs (980425 and 060218) were already observed at distances of 
~ 40 Mpc and ~ 130 Mpc, respectively [26|. The local rate of long GRBs is estimated to be ~ 200 
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Figure 3: Gravitational wave rss amplitude for different accretion rates and Advanced LIGO sensitivity (red 
dashed line) as a function of the gravitational wave frequency, for a — 10° (left panel) and a — 20° (right 
panel). The corresponding outer radius of the accretion disk is indicated in the upper horizontal axis. 



Gpc ^ yr ' (e.g. ^]) and INTEGRAL has detected a significant ratio of faint GRBs which 



are inferred to be local []29[]. As a consequence, this work shows that in the near future the detection 
of precessing disks in the central engines of GRBs should be possible through their gravitational 



emission [3C]. The detection of one event of this class can be used to test the Lense-Thirring effect 
in the strong field limit. 
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